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Abstract 

We consider bipartitions of one-dimensional extended systems whose probability 
distribution functions describe stationary states of stochastic models. We define esti- 
mators of the shared information between the two subsystems. If the correlation length 
is finite, the estimators stay finite for large system sizes. If the correlation length di- 
verges, so do the estimators. The definition of the estimators is inspired by information 
theory. We look at several models and compare the behavior of the estimators in the 
finite-size scaling limit. Analytical and numerical methods as well as Monte Carlo 
simulations are used. We show how the finite-size scaling functions change for various 
phase transitions, including the case where one has conformal invar iance. 
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1 Introduction 



In recent years a lot of attention has been paid to the understanding of the entangle- 
ment properties of ground-state wavefunctions of hermitian Hamiltonians describing one- 
dimensional quantum chain spin systems [U [2] . For a pure state if C = A + B is a bipartition, 
the von Neumann entanglement entropy S v n is defined as: 

S vN (A) = S vN (B) = -Tr(p A \np A ), (1.1) 

where p A = Trs(p), p is the density matrix related to the ground-state wavefunction. If C, 
A and B have L, I and respectively L — I sites, for a large system and large subsystems, S q 
converges to a constant if the system is not critical. If the system is gapless and conformal 
invariant, one gets 

S vN (l,L) ~ 7 lni + C (L»/»l), (1.2) 
and the finite-size scaling behavior 

S vN (l,L) = 7 lnL c + C, L c = Lsm(7il/L)/7i, (1.3) 

where 7 = c/6 for an open system and 7 = 2c/6 for periodic boundary conditions (c is 
the central charge of the Virasoro algebra). C is a nonuniversal constant. The multipartite 
mutual information was also studied [6]- [3] and shown to give more information about the 
system beyond the central charge. 

In the present paper we consider a bipartition of a large one- dimensional interacting clas- 
sical system C and study the shared information between the subsystems A and B due to 
large scale correlations. Such systems can be described by probability distribution functions 
(PDF) of stationary states of stochastic models. They are ket ground-states wavefunctions 
written in a special basis, of non hermitian and in general, nonlocal quantum chain Hamilto- 
nians which give the time evolution of Markov processes. The bra ground-states are trivial. 
The essential difference between the classical and the quantum cases is that in the former 
case the components of the ket wavefunctions are non-negative real numbers which, when 
properly normalized, represent the probabilities to find in the stationary state, various con- 
figurations of the system. 

Wavefunctions with similar properties can also be encountered in the case of some special 
hermitian Hamiltonians. An example is the ground-state wavefunction of the spin 1/2 577(2) 
symmetric Hamiltonian if one uses the valence bond (£77(2) singlets) basis [TJ [8]. 

We present several estimators of shared information which have neat properties. They 
vanish if the subsystems are separated. Their values increase if the shared information 
increases. This happens if the systems are more constrained. It turns out that, in many 
cases, the estimators have a behavior which copies the one of the entanglement entropy. For 
example, if the PDF describes a system with a finite correlation length, the estimators stay 
finite for large systems. The area law [9] is respected. If the system is critical one gets, 
in general, logarithmic corrections to the area law. We are going to discuss in detail these 
corrections. (Some preliminary results were presented in [10].) 

In stationary states, the existence of divergent correlation lengths does not imply confor- 
mal invariance. Therefore one can have a larger variety of critical systems. We are going to 
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show how the estimators behave in several cases. Without entering into details here we just 
mention that if the evolution operator is conformal invariant, the estimators have the form 
(jl.3p with different values for 7. 

Distinct estimators can be sensitive to distinct aspects of shared information. As in the 
case of quantum entanglement, in many cases the evaluation of estimators can be a very 
difficult task. For some estimators however, we can take advantage of the fact that we are 
dealing with stationary states of stochastic processes and use Monte Carlo simulations to 
compute them. 

The estimators of shared information are presented in Sec. 2. The configuration space 
in which the estimators are defined is given by Dyck paths. Dyck paths appear in different 
contexts. They are used in polymer modelling. They can also represent a one-dimensional 
interface. We also recall that if one takes a one-dimensional lattice, and on each lattice site 
one puts a spin 1/2 representation of SU(2), to each SU{2) singlet corresponds a Dyck path. 
This implies that the ground-state of a SU(2) invariant one-dimensional quantum chain can 
be written in a Dyck paths basis. 

Different stochastic models have stationary states described by different PDF of the 
various Dyck paths. In Sec. 3 we present a simple local one-parameter stochastic model (see 
also Appendix A). The stationary states are given by a known model of polymer adsorption 
pTj . Depending on the value of this parameter, one is in a gapped phase (finite correlation 
length) or in a gapless phase which is not conformal invariant. We compute the estimators in 
this model. In Appendix B, we compare the expression of the estimators with those obtained 
considering restricted and unrestricted Motzkin paths instead of Dyck paths. 

In Sec. 4 we consider a different non-local, one-parameter dependent, stochastic model: 
the raise and peel model [121 03] • Depending on the value of the parameter, one is in a 
gapped phase, a conformal invariant phase or in a gapless, non-conformal invariant phase in 
which the critical exponents vary continuously with the parameter. Using analytic methods 
but mostly Monte Carlo simulations, we have computed the estimators and determined their 
expression in the finite-size scaling limit. 

As pointed out in Refs. [T4"| IT5| IT6] in gapless conformal invariant quantum chains, the 
Renyi entropies get different subleading contributions if one consider separately systems with 
/ odd and I even. In Sec. 5 we discuss subleading contributions to several estimators in the 
conformal invariant point of the raise and peel model. 

In Sec. 6 we change the configuration space from Dyck paths to ballot paths. The config- 
uration space is not left-right symmetric anymore. We compute the finite-size scaling limit 
of two estimators in two different models and compare them with the left-right symmetric 
case. 

The estimators for the stationary states of the asymmetric exclusion process (ASEP) on 
a ring are computed in Sec. 7. 

We sum up our results and suggest other applications in Sec. 8. 
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Figure 1: Example of a link pattern for L = 14 and the corresponding Dyck path. In the 
latter there are four contact points and three clusters. The shared information is the largest 
in the left most cluster. The corresponding non-intersecting link pattern (blue lines) is also 
shown. 

2 Estimators of shared information for Dyck paths 
configurations 

We consider an open one- dimensional system with L + 1 sites (L even). A Dyck path 
(restricted solid-on-solid (RSOS) configuration) is defined as follows. We attach to each site 
% non-negative integer heights hi which obey SOS rules: 

h i+1 - ^ = ±1, h = h L = (i = 0,l,...,L-l). (2.1) 

There are 

Z 1 (L)=L\/(L/2)\(L/2 + l)\ (2.2) 

configurations of this kind. If hj = at the site j one has a contact point. Between two 
consecutive contact points one has a cluster. There are four contact points and three clusters 
in Fig. [TJ Contact points and clusters play an important role in shared information. 

In the stochastic processes which we are going to define in the next two sections, the 
Dyck path will be seen as a fluctuating interface between a substrate (hi = (1 + (— l) l+1 )/2; 
i — 0, . . . , L) covered by tiles (tilted squares) and a gas of tiles (not shown in Fig. [TJ. The 
probability distribution functions (PDF) of the various Dyck paths will be given by the 
ground-state wavefunctions of Hamiltonians which give the continuum time evolution of the 
systems. 

There is a known bijection between Dyck paths with L + 1 sites and non-intersecting link 
patterns (SU(2) singlets) for a one-dimensional spin 1/2 system with L sites [17]. In Fig. [1] 
one can see how this bijection is defined in a particular case. This observation opens the 
door for other possible applications of the present work. Spin 1/2 SU(2) symmetric one- 
dimensional quantum chains have their dynamics ruled by Hermitian Hamitonians whose 
ground-state eigenfunctions are in the singlet representation. The XXX chain is an example. 
If one studies the quantum entanglement in this wavefunction, by taking the spin basis, one 
looses the SU (2) symmetry. On the other hand the SU(2) singlet basis is not orthogonal. 
However it might happen that in the basis of SU(2) singlets, the coefficients appearing in the 
wavefunction are non-negative real numbers which through a proper normalisation can be 
interpreted as probabilities [TJ[H]. One can then consider the problem of shared information 
in the wavefunction as if it were the PDF of a stationary state of a stochastic process. 
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Figure 2: A Dyck path {L = 14) with a height hi at the site I, splits into two ballot paths 
a(hi) and b(h{) belonging to the subsystems A and B. 



A bipartition of the system is obtained in the following way. The ensemble of Dyck paths 
(system C with L + 1 sites) is divided into two parts: the sites < % < I (part A) and the 
sites I < j < L (part £>). This implies the splitting of each Dyck path, which at the site I 
has the height hi, into two ballot paths. One RSOS path a(hi) which starts at % — with 
hi = and ends at the site I at the height hi and another one b(hi) which starts at % — I 
with height hi, and ends at % — L with h^ = (see Fig. |2]). 

If we consider the Dyck paths which all have a contact point at % — I (hi = 0), it is 
clear that for these paths, the subsystems A and B are separated. On the other hand, if we 
consider the Dyck paths which reach at i — I a high value of hi, they have around the site 
/ large clusters shared by the two subsystems. This is due to the constraint given by the 
RSOS rules (12. ip . The estimators we are going to define take into account these observations. 
They all measure in different ways the amount of information that can be obtained about 
the ballot paths in B if one observes ballot paths in A. 

We denote by P(a(hi),b(hi)) the probability to have a given Dyck path in C formed by 
the ballot paths a(h{) and b(h{) in A, respectively B. We consider the marginals 

P(a(h l )) = Y,PHhi),b(h l )), (2.3) 

b 

and P(b(hi)). The probability to have a height hi = h at the site / is 

Pt(h, L) = £ P{a(tH)) = £ P{b{h)). (2.4) 

a b 

The proposed estimators are defined as follows: 
(I) Mutual Information, 

/((,£)= £ f(.(tMM) sSi < 2 ' 5 > 
i,«,„(M,i(M P(a(h,))P(b(h,)) 



This is the standard estimator of shared information [18] in information theory. 
(II) Interdependence/ 

This is the Shannon entropy in which the observable is the height at the site I: 

H h (l, L) = - £ Pi(h, L) InPiih, L). (2.6) 

h 
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This estimator is new. It imitates the expression of the von Neumann entropy in which one 
uses the Schmidt decomposition. One can also define in a natural way the Renyi interde- 
pendences: 

R n (l, L) = 1/(1 - n) In £ P t {h, L) n , n = 2, 3, . . . . (2.7) 

h 

(III) Valence Bond Entanglement Entropy 
Is the average height at the site I: 

h(l,L) = Y / hPi(h,L). (2.8) 

h 

This estimator was introduced independently by Chhajiany et al [7] and Alet et al [8] in the 
context of SU{2) symmetric spin 1/2 quantum chains, and further studied by Jacobsen and 
Saleur [TS]. The aim was to measure the average number of link (SU(2) singlets) crossings 
at a given site of the quantum chain. 

We would like to stress that when computing the estimators, the cases of even (odd) 
values of I have to be considered separately at least for small values of I and L. The number 
of Dyck paths which have hi = 1 (/ odd) for example, can be very different than the number 
of paths which have hi = (/ even). This phenomenon affects the values of f l2.4p -( l2T8j) but 
becomes irrelevant for large values of I and L. 

The estimators defined above are sensitive to the probabilities of having large heights 
which implies, due to the RSOS rules, that one has large clusters and therefore large shared 
information. The next two estimators, which are knew, are sensitive to the probabilities to 
have zero heights. 

(IV) Density of Contact Points Estimator 

D(l,L) = -lnP/(0,L) (I even). (2.9) 

This estimator expresses the fact that if the density of contact points p(l,L) = Pi(0,L) is 
small, one has, in general, large clusters. 

(V) Separation Shannon Entropy 

S(l,L) = H(L)-H(l)-H(L-l), (2.10) 

where H(M) = — J2k Pk m Pk ls the Shannon entropy for a system of size M and Pk is the 
probability to have a Dyck path k. The separation Shannon entropy measures the increase 
disorder in C as the effect of putting together the subsystems A and B. All the estimators 
vanish if the subsystems A and B are separated. 

It is easy to see that if all the Dyck paths have the same probabilities equal to l/Zi(L), 

I(l,L) = H h (l,L), (2.11) 

and 

D(l, L) = 8{l, L) = \n{Z x {L)lZ x {l)Z x {L - I)). (2.12) 
To compute (12. lip we use 

Pi(h, L) = d(l, h)d(L - /, h)jZ x (L) (2.13) 
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Figure 3: The stochastic model for polymer adsorption. The tilted square a evaporates with 
a probability q. b is adsorbed with a probability 1 while c is adsorbed with a probability p. 

in (12. 6p . where d(l, h) is the number of ballot paths which start from ho = at i = and 
end at site I with hi = h. If / is even, h also takes even values. It is therefore convenient to 
denote L = 2L', I = 21' and h = 2ti. One has [II]: 

d n> h ') - 2h ' + l 20! 

The relations (12. lip and (12.121) stay valid for any model in which all configurations have the 
same probabilities (see also Section 6). 

In the next two sections we are going to compute and compare the estimators of shared 
information in two very different models. 

3 A stochastic model for polymer adsorption 

We take as the configuration space for the stochastic model, the Dyck paths with L+ 1 sites, 
seen as an interface between a fluid composed of tiles covering the substrate and a gas of 
tiles (Fig. [3]). The model depends on a non-negative parameter u. 

Consider the slope variables = (/ij+i — /ij_i)/2. The evolution of the system in discrete 
time (Monte Carlo steps) is given by the following rules. With a probability Pi = 1/(L — 1) 
each site i of the system is visited. If Si = ±1 the interface (Dyck path) stays unchanged. 
If Si — and hi > h^\ > (local maximum, not part of the substrate), a tile is removed 
(hi — y hi — 2) with a probability equal to q. If Sj = and hi < h^i (local minimum) a tile 
is added (hi — > hi + 2) with a probability equal to p if hi — (the site i is a contact point) 
and a probabiblity to q otherwise. We take p — u, q — 1 (u < 1) and p — 1, q — 1/u — K 
(«>1). 

For u — 1, this is the Rouse model [2D], its extension for other values of u is so obvious 
that probably the model is known. It is a simple exercise to see that in the continuum time 
limit, the Hamiltonian giving the time evolution of the system is non-hermitian except for 
u = 1. 

The ground-state eigenfunction of the Hamiltonian has a very simple expression: the 
components of all the configurations with the same number of contact points are equal and 
are monomials in u~ l . If we take the components of the configurations with no contact point 
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equal to 1, the components of the configurations with m contact points are equal to l/u m . 
(In Appendix A one show the calculation in the case L — 6). The normalization factor of 
the wavefunction is: 

Z K {L) = Y,K m ^\ (3.1) 

M> 

where K = l/u, ip are Dyck paths and m(ip) are the number of contact points of the 
Dyck path. This is the partition function of the well-known model of polymer adsorption 
[2l~j in which K = exp(/3J) is a Boltzmann weight. Let us notice that if u — K — 1, 
all configurations have the same probabilities and the Dyck path can be seen as a one- 
dimensional random walker who starts at the origin and whose position is given by the 
height with time running on the horizontal axis. 

We denote by px{h L) the density of contact points at the site I for a system with L + 1 
sites (L even) and a Boltzmann weight K, and by hx{l,L) the local average height. Their 
averages over all the L + 1 sites are: 

Pk(L) = Wpk(1,L), h K (L) = ]-J2h K (l,L). (3.2) 
L i L i 

Let us sum up some K-dependent properties of the model which can be obtained from 
asymptotic evaluations [2TJ . 



K < 2 



K 



K > 2 



ZM-^^L-**, (3.3) 

PK{L)r.y^L~\ h K {L)~^'\ (3.4) 

Z 2 ~ ^2 L L~^ (3.5) 

p 2 (L) ~ yV 1/2 , UL)-^ 1 ' 2 . (3.6) 



/ x K-2 T , rs K . . 

p K (L) ; r, h K (L) ; r. 3.8 

H V ; 2(A -1)' v ' 2(A-2) v ; 

One notices that in the whole domain < K < 2 one is in the random walker model (K = 1) 
universality class. Surprisingly, the average height has the same values for all K in the whole 
region. Since there are few clusters, one expects the estimators of shared information to have 
large values for large system sizes. 

One has a phase transition at K = 2. At this value of K, the average value of hx{L) 
gets a value which is half the value observed for K < 2 and the density of contact points 
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decreases slower with the size of the system than for < K < 2. Finally, for K > 2, the 
density of clusters is finite which implies that one has many small clusters and therefore one 
expects the estimators to saturate for large system sizes. 

To illustrate how various estimators behave for large system sizes, we take the simplest 
example, taking K = 1. The probabilities of the various Dyck paths are the same and equal 
to \/Zi{L) (see (12. 2p ). We first compute the mutual information which in this case is equal 
to the interdependency (see (I2.1ip -( 12.14p ). It is convenient to take L = 2/7, I = 21' and 
h = 2h' . Using the Stirling approximation one obtains 

P^h',^-—^—, (3.9) 
'L, 
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where z = hi j y L R and 



Ln = HI ~ L = (3.10) 



L 

Notice that Lrw can be written as 

L RW = Lf(x) J f(x) = f(l-x), f(x)~ x (x«l), (3.11) 

with x = l/L and f(x) = x(l — x). We are going to prove that each universality class of 
critical behavior is characterized by a lenght L with a distinct function f(x). From (12.61) 
and (13.91) we get 

/(/, L) = H h (l, L) ~ ~ ln L RW + C, (3.12) 

where 

C = i -- In --- w 0.303007. (3.13) 
'2 2 2 v ; 

7' is the Euler constant. 

The Renyi interdependency ( 12. 7ft has the same asymptotic expression as the interdepen- 
dency: 

Rn(l,L) ~ l/21nL R + C n . (3.14) 
From (13. 9 p we also get the expression of the valence bond entanglement entropy (12.81) : 

h(l,L)~-4=L)gr. (3.15) 



2vr 

Finally, one can compute the separation Shannon entropy, equal in this case to the density 
of contact points estimator, using ( 12 . 1 2[) and the Stirling approximation 

D{1, L) = S(l, L) ~ 1 ln L RW + ln (3.16) 

The local density of contact points being: 

PK«> L )~^JW- (3 - 17) 
V L RW 
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We observe that in the finite-size scaling regime all the estimators depend on the variable 
Lji, which for large values of L' becomes equal to I'. This is similar to the finite-size scaling 
variable Lsin(7r//L)/7r appearing in section 5. 

We first notice that the valence bond entanglement entropy h(l, L) behaves as a power 
for large distances. The other four estimators have the following behavior in the finite-size 
scaling limit: 

E{l,L) ~ lE \nL RW + C E , (3.18) 

where for 1 << / << L, 

E(l,L) ~ lE \nl + C E , (3.19) 

with 7/ = jjj = 1/2 and 7r> = 7s = 3/2. 

An obvious question is: are the divergent part of the estimators of shared information 
universal? This would imply that in the whole domain < K < 2, in the expressions (13. 18ft - 
(I3.19P only the constants C E change. This was checked using Monte Carlo simulations for 
the interdepedency as well as for R 2 (l,L). In the case of the valence bond entanglement 
entropy h(l,L), this implies that (I3.15P is independent of A. 

~ 1/2 

If one assumes for general A the finite-size scaling form h(l, L) = ckL^ w than from the 
A-independence of hx in f l3.4[) we get that c# is independent of A, and hence universality. 
The same argument can be used in the case of D(l, L) using (13.17[) and the A-independence 
of the power law fall off with the size of the system of p(K, L) given by ( 13.4ft . 

We show now that for large systems, the separation Shannon entropy is also universal. 
The Shannon entropy for a system of size L is 

K m I K m \ 

^ m A K [L) \Z K {L)J 

where is the probability of a Dyck path if), D(L, m) is the number of Dyck paths with m 
contact points, and 

Z K (L)=J2K m D(L,m). (3.21) 

m 

From (I3.20p and (I3.2ip one obtains: 

H K {L) = In Z K {L) - Lp K (L) In A. (3.22) 
We now use (13.31) and (13 ,4p to get the Shannon entropy for a system of size L 

H K (L) = Lln2--lnL-^^lnA + ln I - 2K - J-) . (3.23) 
Ky 1 2 2- A ^(2- A) 2 V tv J K J 

The 7eparation Shannon entropy (12.101) is: 



S K (l,L) = H K (L)-H K (l)-H K (L-l) 

= l>»^ + |z#^->»(p^l/rl. (3-24) 
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where Lrw = 1(1 — |-). We can conclude that the leading term is universal. 

There is a supplementary argument for the universality of the coefficient 75- = 3/2. If 
instead of Dyck paths one takes, with equal probabilities, Motzkin paths (see Appendix B 
for their definition), the asymptotic expression of the partition function is: 

Z " (L) = \/& (3 ' 25) 

with the same leading behavior as in f)3.24j) : 

S K (l,L)~^\nL RW + ±\n(^p). (3.26) 

This is to be expected since one has again a random walker who, at each step, can not only 
move up or down but can also stay in place. This corresponds, in the continuum, only to a 
change of the diffusion constant. More interestingly, in Appendix B we consider restricted 
Motzkin paths. They correspond to SU(3) link patterns [22]. Because of the constraint one 
expects more shared information and a larger value for 75. This is indeed the case. One 
finds 7s = 4. 

We did not compute the mutual information for values of K ^ 1 but we assume that one 
gets 7/ = 1/2 for the whole interval < K < 2. 

To sum up, in the whole critical domain all five estimators of shared information have a 
universal, albeit different, large lattice sizes behavior. We did not study properly the K = 2 
case. Heuristic arguments suggest that the shared information is smaller than in the case 
of the random walker universality class (K < 2). Comparing h,2(L) given by ( 13. 6 j) with 
Hk{L) given by ( 13 AL for K < 2, we conclude that the valence bond entanglement entropy 
for K = 2 is half of the K < 2 value. From the expression ( 13. 6 h of the average density of 
contact points one expects a value 7/5 = 1/2 as compared with 7^ = 3/2 for K < 2. 

We discuss now the K > 2 domain in which the correlation length is finite. In this 
domain one has a finite density of small clusters and a finite averaged height (see ( I3.8P ). 
This implies that all the estimators have a finite value for large lattice sizes. This can be 
seen directly for the valence bond entanglement entropy and for the density of contact points 
estimator by using (I3.8p . and the separation Shannon entropy by using ( I2.10p . (13. 7p . (I3.8P 
and ( I3.22p . The finiteness of the remaining two estimators comes from the exponential fall 
off of the probabilities for large size systems. 



4 Shared information in the Raise and Peel model 

This model has been studied in detail in several papers ( |12[ [23] . for a review see [13]). We 
are going to mention here only those facts which are relevant for the study of the shared 
information. 

Like in the previous section, we consider the interface separating a film of tiles deposited 
on the substrate, from a rarefied gas of tiles. The model depends on a non-negative parameter 
u and the evolution of the system in discrete time (Monte Carlo steps) is given by the 
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Figure 4: Example of a configuration of the interface in the raise peel model for L = 18. 
Depending on the position where the tilted tiles reach the interface, several distinct processes 
occur (see the text). 



following rules. With a probability Pj = 1/(L — 1) a tile from the gas hits the site i = 
1,...,L-1. 

Depending of the slope Sj = (hi+i — hi-\)/2 at the site i, the following processes can 
occur: 

1) Si — and faj > fai-i- The tile hits a peak and is reflected. 

2) Sj = and hi < hi-\. The tile hits a local minimum and is adsorbed (hi — > hi + 2) with 
a probability p and is reflected with a probability 1 — p. 

3) Si = 1. The tile is reflected after triggering, with probability q, the desorption (/ij — > hj—2) 
of a layer of b — 1 tiles from the segment {j=i+l,...,i + £>— 1} where hj > hi = hi + b (see 
Fig. |4]). With a probability 1 — q nothing happens. 

4) Si = —1. With probability q, the tile is reflected after triggering the desorption {hj — > 
hj — 2) of a layer of b — 1 tiles belonging to the segment {j — i — b + 1, . . . , i — 1} where 
hj > hi = hi_ b . With a probability 1 — q nothing happens. 

We take p = u, q = 1 for u < 1 and p — 1, q = 1/u for u > 1. In Fig. 4 we illustrate the 
possible processes in the case of a Dyck path with L = 18. The tilted tile b corresponds to 
process 1), tile c to process 2), tile a to process 3) and tile d to process 4). 

The phase diagram of the model was determined from the study of the finite-size scaling 
limit eigenspectrum of the Hamiltonian which gives the time evolution of the system, and 
other methods. For < u < 1, the system is massive. For u — 1, the system is conformal 
invariant with an effective central charge of the Virasoro algebra c = 0. For this case only, the 
stochastic process is described by a spin 1/2 XXZ quantum chain with U q (sl(2)) symmetry 
with q = exp(i7r/3). The link pattern corresponds to U q (sl{2)) singlets [T7]. For u > 1, 
the system is scale invariant but not conformal invariant. The dynamic critical exponent z 
decreases continuously from the value z = 1 when u — 1 to very small values when u becomes 
large. What makes this model special is the fact that it contains at u = 1 a stochastic model 
which is conformal invariant and we can therefore compare the behavior of the estimators 
with those studied in the problem of quantum entanglement. The dynamical properties 
of the model are reflected in the properties of the stationary state. For < u < 1, the 
densities of clusters are finite for large system sizes (like in the polymer adsorption model 
for u < 1/2). For all values of u > 1 the density of clusters vanishes with a power which 
varies continuously with u. This was not the case in the polymer adsorption model where 
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128 



Figure 5: Typical configurations in the stationary states for L + 1 = 129 sites and three 
values of u. 



for all values of u > 1/2, one is in the random walker phase. Some typical configurations for 
three values of u are shown in Fig. [5] for a lattice of L = 128 (129 sites). One notices many 
clusters for u = 0.4 (in the massive phase), fewer clusters for u = 1 and a single cluster for 
u = 2.5. One finds fewer and fewer clusters when u increases. 

We discuss the shared information in the stationary states of the model. In the domain 
< u < 1, the density of clusters is finite for large values of L and so are the average values 
of the heights. As discussed in the previous section, this implies that all the estimators stay 
finite for large L. The situation is different when the system is critical. 

We first discuss the u = 1 case. We will show that all the estimators of shared information 
have the following expression in the finite-size scaling limit: 

E(l,L) ~ lE i n l c + Ce , (4.1) 

with 

L c = Lsin(7r//L)/7r, (4.2) 

from which follows 

E(l,L) ~ lE \nl + C E , l«l«L. (4.3) 
These expressions are similar to what was found for the von Neumann entanglement entropy 

(OD. 

The mutual information (12. 5p was obtained using the exact PDF of the stationary states 
up to L = 26 and assuming the finite-size scaling behavior ( 14. ip . The data are compatible 
with ( 14. ip if 7/ = 0.07 and Cj = 0.65. This is a useful estimate rather than a result. 

In order to proceed, using Monte-Carlo simulations, we have determined the probabilities 
Pi(h, L) of having a height h at the site / for a lattice of size L and stored the information. 
This is possible since even for L = 40000 the probability to obtain a height larger than 70 is 
negligible. These data have been used for several purposes. 

We first studied the interdependency Hh(l,L). To determine 7# we have taken large 
lattices (up to L — 40000) and 1 << I « L. In Fig. [6] we show the plot Hh(l,L) as 
a function of In / for various lattice sizes. As expected we have data collapse. The larger 
the values of L, the larger is the domain of data collapse. We have fitted the L = 40000 
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0.8 



0.7 



L=40000 
L=30000 
L=20000 
L= 10000 
L= 4000 



In (I) 



Figure 6: The interdependency Hh(l,L) as a function of In / for various lattice sizes L = 
4000,10000,20000,30000,40000. Also shown the two domains where two linear fits were 
made in order to determine 7# and Ch- 



data assuming straight lines in two regions shown by the two squares of Fig. |6j One obtains 
H h = 0.055 In Z+0.62 for the low values of In laMH h (l,L) = 0.051 In Z+0.65 for the larger In/ 
domain. We take the result of the second fit for our estimates of 7# and Ch since probably, 
in the first domain the condition 1 << I was not satisfied. Maybe not surprisingly, the values 
for 7# and Ch obtained for the interdependency are close to those obtained for the mutual 
information 7/ and Cj. As shown in Sec. 3, if the probabilities for all configurations are equal 
(this is not the case here), the mutual information and the interdependency are equal. In 
order to find the expression of the interdependency in the finite-scaling limit, we have taken 
a large lattice [L = 40000) and computed the quantity (H h (L/2, L) — H h (l, L)) as a function 
of ln(sin(7r//L)). If the finite-size scaling expression (14.1 j) is correct, one should obtain a 
straight line with a slope equal to jjj. In Fig. [7] we show the results of the Monte Carlo 
simulation and one can see that in the finite-size scaling limit, Hh{l/L) has the expression 
(14. ip for u — 1, which is not the case for u — 4. 

The Renyi interdependency R n (l,L) were also measured. We have found 7r 2 = 0.05, 
Cr 2 = 0.4 and 7# 3 = 0.04, Cr 3 = 0.4. We have to stress that a precise determination 
of the values of 7# n for the Renyi interdependencies using Monte Carlo simulations is not 
an easy task. The values of I, for which one can assume that the finite-size scaling limit 
is valid, are large. This can be seen from the non-negligible differences between the Renyi 
interdependencies computed for I odd and I even due to subleading terms in (14. ip . In Sec. 5 
we discuss in detail these subleading contributions. 

The valence bond entanglement entropy h(l, L) was also measured. We have found 7^ = 
0.277 and Ch = 0.75. The determination of 7^ and the check of the finite-size scaling function 
(14. ip are much more precise in the case of h(l, L) as shown in Figure 2 of [10]. Non-leading 
effects play a smaller role. Considering periodic boundary conditions and the boundary 
Coulomb gas formalism, Jacobsen and Saleur [12] have computed 7^ in a related model. 
Their value divided by 2 (we have an open system while in their case the chain is periodic) 
is 7 = V3/2n w 0.275, very close to the measured value 7^. This observation is not trivial 
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Figure 7: H h (L/2,L) 
(L = 40000). 



In[sin(%l/L)] 

— Hh(l,L) as a function of ln(sin(7r//L)) for u = 



1 and u = 4 



because it proves that valence bond entanglement entropy is a bona fide estimator. In the 
periodic system the information is shared at the two ends of the subsystems. 

Moreover, we have determined the second cumulant k 2 of the probability distribution 
Pi(h, L), it has the same finite-size scaling behavior as the estimators: 



We found k 2 = 0.19 and b 2 = 0.25 is in excellent agreement with half the value obtained by 
Jacobsen and Saleur [19]: (27iV3 - 9)/vr 2 w 0.190767. The fact that k 2 in our case (open 
system) is half the values obtained in [19], which corresponds to a periodic system, indicates 
that k 2 can also be an estimator for the shared information between subsystems. 

The density of contact points estimator was computed using an almost rigorously derived 
expression [23] for the density of contact points: 




k 2 (1,L) ~ (3 2 \nL c + b 2 . 



(4.4) 




p(l,L) 



a 



(4.5) 




where 




(4.6) 



Using ( 12. 9 p we have: 



D(l, L) In p(l, L) = - hiL c + 0.28349. 

3 



(4.7) 
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The separation Shannon entropy was estimated using the same data as for the mutual 
information. We found 75 = 0.4 and C$ = 0.7. 

We notice that the values of 7^ for the density of contact points estimator and 75 for the 
separation Shannon entropy are close to each other but very different of those determined 
for the mutual information and the interdependency. One has to keep in mind that D(l,L) 
and S(l, L) coincide if the configurations have the same probabilities. 

To sum up, in the conformal invariant point of the model, all estimators have the expres- 
sion (14.11) with several values for 7 (see Tabled]). 





u = 1 


u = 1 


u = 4 


u = 4 




lE 


C E 


lE 


C E 


Mutual information 


0.07 


0.65 






Interdependency 


0.050 


0.67 


0.09 


0.91 


Renyi (n — 2) 


0.05 


0.39 


0.06 


0.09 


Valence bond ent. 


0.277 


0.75 


0.63 


1.37 


Dens. Contact points 


0.333 


0.284 


0.73 


0.71 


Separation Shannon 


0.4 


0.7 







Table 1: The values of 7^ and C E in QH) and fT4T2|) 



We have studied the estimators of shared information for several values of u > 1. We 
are presenting the results for u = 4 and just mention what happens at other values of u. 
For u = 4 the system stays critical, the dynamic critical exponent is z ~ 0.3 and conformal 
invariance is lost (z ^ 1). 

From the Monte Carlo simulations, the following picture has emerged. Equation (14. 3 j) 
stays valid but in the finite-size scaling limit one has, instead of (14. 2ft . 

E(l,L)~ lE lnL u=A + C E , (4.8) 

where the L u=i = Lf u=i {l/L) function has replaced the function L c = L u=1 . That one has 
a new finite-size scaling function can be seen, for example, from Fig. [7] and from Figure 2 of 
|10j . The fact that there is only one finite-size scaling function for all the estimators can be 
seen, for example, from Fig. |HJ Taking different values oil/L, we have plotted the differences 
between the values of two estimators (H(l, L) and h(l, L)) at / = L/2 and I. If there is only 
one scaling function one should obtain a straight line with a slope given by the ratio ^h/lu- 
One can see that for large values of L, this is the case. 

The values of 7^ obtained for u = 4 are given in Table 1. We notice that they are all 
larger when compared with the values at u = 1. The shared information is larger. We notice 
that the values of Ce also increase. 

We have explored also different values of u. The larger the value of u, the larger are the 
values of 7^ and Ce- The shape of the finite-size scaling function L u changes continuously 
with u. 



16 




L=40000 
L=30000 
L= 15000 
L= 10000 
L=4000 



0.05 



0.1 0.15 0.2 
H h (L/2,L)-H h (l,L) 



0.25 



0.3 



Figure 8: (h(L/2,L) — h(l,L)) as a function of (Hh(L/2, L) — Hh(l,L)). The data were 
obtained from lattices sizes L = 4000, 10000, 15000, 30000 and 40000. 



5 Parity effects and subleading contributions to esti- 
mators in the raise and peel model at the conformal 
invariant point. 

As noticed in Sec. 4, for small lattices there are important differences between the Renyi 
interdependencies computed for I odd and I even at u — 1. A similar phenomenon was 
observed for the Renyi entropies but not for the von Neumann entanglement entropy, in the 
study of quantum chains [Til E2 EE]- Moreover, in [TJ] it was conjectured that for the XXZ 
model with an anisotropy — 1 < A < 1 and open boundaries, in the finite-size scaling limit, 
one has: 

6R n (l,L) = f n (l/L)/L^ n (5.1) 

where 

5R n {l,L) = R n (l,L) - R n (l + 1,L) (/even) (5.2) 
i? n (/,L) = T ^- In Tr(p^) (5.3) 

and 

K = 7r/(2arccos(A)). (5.4) 

Pa is defined in (11. ip and L c in (14. ip . In ( 15 .ip / n (0) is finite. This conjecture was checked 
analytically at the decoupling point (A = 0) and numerically for other anisotropics. For 
periodic boundary conditions, K has to be replaced by 2K in (15. ip . 



17 



We have decided to look to the odd-even finite-size corrections for some of the estimators 
discussed in Sec 4 and see if one gets a behavior similar to (15.11) . We have considered the 
quantities 

6E{l,L)=E(l,L)-E{l+l,L) (/even) (5.5) 

where E(l,L) is an estimator. 

In the case of the interdependency, the data shown in Fig. suggest an expression similar 
to fl£H): 

5H ~ ci/Z* 1 (5.6) 

with c\ ~ 1.3 and x\ ~ 0.58. The value of the exponent X\, as expected, does not match 
the exponent K = 3/2 (n = 1, A = 1/2 (12] ) in (15. ip . The bases in the quantum chain 
and in the raise and peel model are completely different and there is no reason to have any 
connection between the exponents. 
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Figure 9: Odd-even corrections to the interdependency (see (15.51) and (15. 6p . We have used 
lattices of sizes L = 4000, 20000 and 40000. 



A new phenomenon, not seen is the quantum chains, appears if one considers the Renyi 
interdependencies ( 12. TP . We have measured 6R2U, L) and 8R${1, L) using lattices up to 40000 
sites. The data shown in Figs. [10 and HH suggest the following expressions for the odd-even 
corrections: 

5R n ~ l/L x c n F n (\n(L c )) (5.7) 

Where Xi ~ 0.45, X3 ~ 0.32 and F n (y) are oscillating functions. This observation came as 
a surprise. It is important to stress that this is not an effect of the finite-size scaling limit 
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(l/L finite). Oscillations are seen also in the limit L very large and I finite. Notice that like 
in fl5.ll) we have X\ > x 2 > x 3 . 

We have no explanation for this oscillatory behavior and it would be interesting to know 
if this is a pathology of the raise and peel model at the conformal point or the phenomenon 
is more general. We have looked at the odd-even corrections at higher values of u. The cor- 
rections become small and it is hard to make statements based on Monte Carlo simulations. 




Figure 10: The odd-even corrections for the R2(l,L) Renyi interdependency multiplied by 
L° 45 , a a function of lnL c for different lattice sizes (L = 4000, 20000 and 40000). 



We have also looked at the odd-even corrections for the valence bond entanglement 
entropy (12.81) and found (see Fig. [T2l 

5h(l,L) = c h /L x c h (5.8) 

with Ch ~ 0.50 and Xh ~ 0.99. No oscillatorry behavior was seen. 

Odd-even effects are also seen for the estimators of shared information in the case of the 
polymer adsorption model but we are not going to discuss them here. 
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Figure 11: The odd-even corrections for the R^(l,L) Renyi interdependency multiplied by 
L°- 32 , a a function of lnL c for different lattice sizes (L = 4000, 20000 and 40000). 



6 Shared information in non left-right symmetric sta- 
tionary states. 

The two stochastic models considered in the last sections were left-right symmetric. As a 
consequence, the expressions of the estimators in the finite-size scaling limit were of the form: 

E(l,L) ~ lE \n(Lf(x)) + C E (6.1) 

where f(x) = f(l — x) and fix) ~ x, for small values of x. In this Section we consider 
a stochastic process with a source at the left-end of the system in order to see how the 
asymmetry in the system is reflected in the estimators of shared information. We will 
consider two examples only. 

Instead Dyck paths, we consider ballot paths. Those are paths respecting the RSOS 
rules (12. ip with a modification: Jil is fixed like for Dyck paths (hi = 0) but ho is free 
(ho = 0,2, . . . , L). See Fig. 9 for an example. There are L\/((L/2)\) 2 configurations of this 
kind. 

Assume that each of these configurations has equal probability. One can easily compute 
the finite-size scaling limit of the density of contact points estimator (equal to the separation 
Shannon entropy), one obtains: 



D(l,L) = S(l,L) ~ - In 1(1 



z><z> 1/3 
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Figure 12: Odd-even corrections to the valence bond entanglement entropy. Lattices of sizes 
L = 4000, 20000 and 40000 were used. 
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Figure 13: A ballot path for L = 14. One has one contact point and one cluster of size 
14. 

= §lnW-l»(l) + iln(f). (6.2) 

In deriving this result we have taken ballot paths between the sites and I (hi = 0) and Dyck 
paths between the sites I and L (fi£ = 0). One can compare the expression ( 16. 2ft with (13.161) 
(the same estimators for Dyck paths). One notices that the value of 7^ = 75 is the same but 
instead of having a constant like in (I3.12p . one has an asymmetric (under the interchange 
l/L -H- 1 — l/L) function of l/L. If l/L approaches the value 1, one recovers the expression 
(I3.13P as one should since one is far away from the "source" at the site i — 0. The limit 
l/L — > is singular since for L large there are no contact points for I close to the origin. We 
are going to find a similar situation in our next example. 

The raise and peel model with a wall (RPMW) is a simple extension [23] of the raise and 
peel model described in Sec. 5 . We consider this extension for the case u = 1 only. This is 



21 




oi 6 n ;,= m oi e ii /.= ii 



Figure 14: The adsorption of a half-tile at the first site and desorption of a layer touching 
the boundary for the ballot path shown in Fig. 9. 
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Figure 15: Typical configuration in the stationary state for the RPMW with a boundary rate a = 1. 
The system has a size L = 128. hi is the height at site i. The contact points are concentrated at 
the right end. 



the case where one has conformal invariance. 

The stochastic model is denned by the four rules given in Sec. 5 to which we add the 
following rules: with a probability Pj = 1/(L + a — 1) a tile from the gas hits the site i, 
(i — 1, . . . , L — 1). The changes of the interface produced by the hits are the same as in the 
raise and peel model. With a probability Pq = a/(L + a — 1), a half-tile hits the site 0. The 
boundary rate is equal to a as opposed to the bulk rates which are equal to 1. The "source" 
rate a stays as a free parameter in the problem. If the slope s — hi — h is equal to 1, the 
half-tile is adsorbed (see Fig. [6]). If Sj = — 1, the half-tile is reflected (see Fig. 10). 

The properties of the stationary states of the raise and peel model with a wall have been 
described in detail in [23]. We are interested here in one result only. It was found that in 
the finite-size scaling limit, for any boundary rate a, the density of contact points has the 
expression: 

(\ _ cosf^U 1/3 

* L »~ 8757 ) ■ (M) 

From which we derive the following expression for the density of contact points estimator: 



D(l,L) ~ -InLc + ln 
3 



0.8757 f 1 - cos A 



1/3' 



(6.4) 
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One notices that for l/L close to 1, one is back at the result (14. 7p obtained in the raise and 
peel model and similar to Eq. (16. 2p . the limit l/L — > is singular. We didn't look at other 
estimators. 



7 Shared information in the asymmetric exclusion pro- 
cess with periodic boundary conditions. 

To illustrate how the estimators defined in Section 2 can also be used for configurations 
which are not Dyck paths, we give a simple example: the asymmetric exclusion process on 
a ring. 

We do not describe here the very well known asymmetric exclusion process (see [21] for a 
review). We consider directly the stationary states. For a ring with L sites and N particles, 
in the stationary state all 

i\n 

Z(N,L) = „ — (7.1) 

v ' ; N\(L-N)\ v ; 

configurations have the same probability. We are interested in the shared information be- 
tween a segment with I sites and the remaining segment with L — I sites on the ring. The 
average local density is r = N/L. The shared information between the two segments comes 
from the fact that if we measure Q particles in the segment /, one has to have N-Q particles 
in the segment L — I. 

The mutual information, equal to the interdependency (see ( 12. 12ft ). has the expression: 

1(1, L) = H q (l, L) = -Y, Pi(Q, L) InPtiQ, L), (7.2) 

Q 

where 

Pi(Q, L) = Z(Q, l)Z(N -Q,L- l)/Z(N, L). (7.3) 

Denoting Q = rl + q and taking \q\/L << 1, in the finite-size scaling limit, Pi(q, L) is given 
by a Gauss probability distribution: 

%L)~_L=exp(-^l), (7.4) 

where 



a = yV(l -r)L RW , (7.5) 
and Lrw is given by (13. lip . From ( 17.2ft and (17.41) one obtains. 

1(1, L) = H q (l, L) ~ ~ In L RW + ~ ln(2vrr(l ~r)) + \. (7.6) 
The equivalent of the valence bond entanglement entropy is: 

q(l, L) = Y, L)~2a = 2^(1 - r)L^. (7.7) 



23 



Finally, the equivalent of the density of contact points estimator, equal to the separation 
Shannon entropy is: 



We notice that the expressions of the estimators are very similar to those obtained in 
the case of the model of the polymer adsorption in the domain < K < 2. One difference 
is that in the present case, four of the estimators have the same leading expression and that 
the common coefficient 7 # is smaller (one has to take half the value obtained in the present 
case since we have periodic boundary conditions). Moreover, half of the coefficient in front 
of L)(w in ( I7.8P is smaller than the coefficient of L 1 ^ in (I3.15p . This implies that the shared 
information is smaller than in the polymer adsorption model. 

~ 1 /2 

Another difference is that the coefficient in front of L^ w in (17. 7p depends on the density 
of particles r. The 7^'s coefficients (equal to 1/2) are not. The r-dependence being hidden 
in the constant term. 

8 Conclusions. 

Our aim was to see if for stationary states of one-dimensional Markov processes with many 
degrees of freedom, one can define estimators of shared information between two subsystems. 
We wanted to make a connection between information theory and non-equilibrium statis- 
tical physics paralleling the connection between quantum information theory and extended 
quantum systems. We found several useful estimators. The Markov process gives the rules 
how to perform Monte Carlo simulations. These simulations on large system sizes were used 
to compute some estimators. This method can not be used in the quantum case. 

Although we consider classical systems, it turns out that the estimators have properties 
similar to their counterparts in the quantum mechanical case. 

The estimators are defined in Sec. 2. We have studied their properties in several models 
described in Sections 3, 4, 6 and Appendix B. If the subsystems have lengths / and L — I, 
all the estimators have the following behavior in the finite-size scaling limit (1 << l,L, l/L 
fixed) : 

a) They vanish if the subsystems are separated (for any / and L). 

b) If the correlation length is finite, they stay finite for large values of L (area law). 

c) If the system is critical, for each universality class one can define a characteristic length 



where x = l/L and f(x) depends on the universality class. 

d) At criticality one obtains corrections to the area law. They all depend on L only and 
can increase logarithmically or as a power of L. 

e) If one has logarithmic corrections, the estimators E(l, L) have the following expression: 







E{l,L) ~ 7E In L + C E , 




where 7 # depends on the estimator but is universal, Ce is not universal. 
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f) If one has power-law corrections, the estimators have the following expression: 

E(l,L)~5 E (L) m + D E , (8.3) 

where the exponent m < 1 and the constant S E are universal and D E is not universal. 

g) If one compares the coefficients ■y E and 5 E belonging to two universality classes, they 
are all larger (smaller). 

In our study we have found only one estimator (the valence bond entanglement entropy 
( 12. 8p ) for which one gets power-law corrections with m = 1/2. This was seen in the models 
discussed in Sec. 3 and 6 but not in the model of Sec. 5. All the other estimators give 
logarithmic corrections. We don't yet have an explanation of this observation. 

We do have some understanding of the physical ingredients which enter the expression 
of the estimators. 

The simplest case is when we have conformal invariance (u — 1 in the raise and peel 
model of Sec. 4). The expression of L = Lq ( 14. 31) is the same as in the quantum case and 
is well understood pQ. The value of 7^ (see Table 1) can be inferred from the calculation of 
Jacobsen and Saleur [19] for the periodic case, is related to the critical exponent of the 
local density of contact points which is known exactly. The values of 7/, 7^ and 75 have 
still to be derived analytically. Subleading contributions to ( 18. 2 j) observed if one considers 
even or odd sites are discussed in Sec. 5. In the case of the Renyi interdependencies, a new, 
unexpected phenomenon is seen: the differences between the Renyi interdependencies for 
odd and even sites oscillate. 

We do not have a clue how to derive the values of j E s or the functions L u for the whole 
domain u > 1 of the raise and peel model. 

It looks as if the models of Sec. 3, 6 and Appendix B have one common L with f(x) = 
x(l — x). These models share a common feature: the heights respectively the hopping 
particles have Gauss fluctuations. 

We have also studied simple examples of left-right asymmetric partitions and found the 
finite-size scaling functions in these cases. 

The estimators defined here for bipartitions can be extended in an obvious way to mul- 
tipartitions. The simplest case is the separation (contact) points of tripartitions, which is 
related to a not connected two-point correlation function. We did not touch this topic in 
this paper. 

We hope to see the application of the methods introduced here to other stochastic pro- 
cesses. Probably new aspects of the shared information theory will unravel. 
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|4> |5> 
Figure 16: The five Dyck path configurations for L = 6 



A Polymer adsorption for seven sites 



In this appendix we derive the Hamiltonian and the expression of the stationary state for 
the model described in Section 4 in the case L = 6. 

There are 5 Dyck paths configurations in this case shown in Fig. [161 The time evolution 
of the probabilities P a {t) to find the system in the state \a) (a = 1, 2, . . . , 5) at the time t is 
given by the 5x5 matrix H a ^. Its non-diagonal matrix elements are given by the rates b — > a 
with the sign changed {H a ^ < 0) and the diagonal elements are fixed from the conditions: 
Y, a Ha,b = 0. 

Denoting 



\P(t)) = J2 p a(t)W), Pa = lim P a (t), |0> = £P» 



a=l 



0=1 



we get 



d/dt\P(t)) = -H\P(t)), H\0) = 0. 
The Hamiltonian with the transitions rates fixed by the rules given in Section 4 is: 



(a.i; 



(A.2) 



H 
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(A.3) 



The stationary state follows: 

|0} = |l) + |2) + u- 1 (|3) + |4}) + u- 2 |5). 



(A.4) 



B Restricted and unrestricted Motzkin paths. 

A Motzkin path is defined as follows. Consider the Cartesian plane ZxZ with the sites 
on the x-axis and the heights on the y-axis. The lattice paths start at (x, y) = (0, 0) and 
one can use steps (U', L', D'), where U' = (1, 1) is an up-step, L' — (1, 0) is a level-step and 
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D' = (1,-1) is a down-step. The steps end at (x, y) = (L, 0) without going below the x-axis. 
There are 

L/2 D 

Zm{L) = £ k\(k + 1)\(L - 2k)\ (B ' 1} 

configurations of this kind. For large values of L one obtains (13.251) . 

If one takes the Zm{L) Motzkin paths with equal probabilities, they describe a random 
walker who starts and returns to the origin after L steps moving with a diffusion constant D 
differing from the Dyck paths cases. This explains (13.260 . The leading finite-size behavior 
of the estimators is independent of D. 

The restricted Motzkin paths are defined for L = 3n and are subjected to the constraint 
that for any < x < L one has Nu'(x) > Nl>(x) > No'{x). This means that after x steps, 
the number of up-steps Njj/ has to be larger or equal than the number of level-steps Ny 
which are larger or equal the number of down-steps Nr>>- 

The number of restricted Motzkin paths is equal to [22] : 

2L! 

Z RM {L) = — ( — —-T-, — — -t. (B.2) 
n\[n + ly\n + 2)! 

There are obviously fewer restricted Motzkin paths than unrestricted ones. For example if 
L = 3 one has only one restricted path as opposed to four unrestricted ones. The number of 
restricted Motzkin paths is equal to the number of independent s/(3) singlet configurations 
one gets in a one-dimensional lattice with L+l sites, if one puts on each site the 3-dimensional 
fundamental representation of sl(3). 

Due to the fact that restricted Motzkin configurations are more constrained than the 
Dyck paths, one expects in the case of a bipartition, larger shared information between the 
two subsystems. This is indeed the case. Using the (12.121) with Zi(L) replaced by Zrm one 
obtains the following finite-size scaling expression for the separation entropy (equal, in this 
case, to the density of separation points estimator): 

S(p,n) = 41n[p(l -p/n)]+ln(7r/V3), (B.3) 

where we have taken L = 3n and / = 3p the number of sites in the bipartition. Comparing 
(IB.3j) to (I3.16P we notice that the finite-size scaling function is unchanged but the value of 
75 jumped from 3/2 for Dyck paths to 4 in the case of restricted Motzkin paths. 
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